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Abstract. We present a novel determination of the astrophysical uncertainties associated to 
the secondary antiproton flux originating from cosmic-ray spallation on the interstellar gas. 
We select a set of propagation models compatible with the recent B/C data from PAMELA, 
and find those providing minimal and maximal antiproton fluxes in different energy ranges. 
We use this result to determine the most conservative bounds on relevant Dark Matter (DM) 
annihilation channels: We find that the recent claim of a DM interpretation of a gamma-ray 
excess in the Galactic Center region cannot be ruled out by current antiproton data. 

Einally, we discuss the impact of the recently released preliminary data from AMS-02. 
In particular, we provide a reference model compatible with proton, helium and B/C spectra 
from this experiment. Remarkably, the main propagation parameters of this model are in 
agreement with the best fit presented in our earlier statistical analyses. We also show that 
the antiproton-to-proton ratio does not exhibit any significant anomaly at high energy with 
respect to our predictions. 
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1 Introduction 

The quest for some anomaly in the cosmic-ray (CR) antiproton spectrum has captured the 
interest of the astro-particle community for more than two decades and is presently one of 
the main targets of the AMS-02 observatory [1]. 

The rationale of this effort is the following: If the astrophysical sources, mainly super¬ 
nova remnants (SNRs), do not inject primary p in a relevant amount ” and provided that 
the secondary antiproton flux due to CR spallation on the diffuse interstellar gas can be ac¬ 
curately computed - then the measured CR antiproton spectrum becomes a valuable probe 
of new physics (see, e.g., [2] and references therein). In particular, a fascinating scenario is 
that p are copiously yielded by WIMP dark matter (DM) particles annihilating in the dark 
halo of the Galaxy [3]. 

In this context, an important role is played by secondary-over-primary ratios of CR 
nuclear species: Most noticeably the boron-to-carbon (B/C) ratio is one of the most useful 
tracers of Galactic CR propagation since boron is expected to be not produced in stars and 
the cross sections for B production from its main primaries (C and O) are better known. This 
channel offers a rather robust constraint on CR models and allows a precise determination of 
the secondary!? flux produced in the Inter-Stellar Medium (ISM) (see, e.g., [4-8]) Moreover, 
it is crucial to constrain the amount of antiprotons produced and then accelerated in SNRs [9, 
10 ]. 
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In this paper we make use of the recent measurements of the B/C ratio, together with 
the proton, helium and light nuclei absolute spectra provided by the PAMELA collabo¬ 
ration [11, 12] to constrain the propagation parameters, and to accurately compute the 
secondary antiproton spectrum in the 0.1 — 10^ GeV energy range. 

Our main goal is to achieve a comprehensive evaluation of the relevant uncertainties 
affecting this computation, extending the work presented in [13]. We focus our attention 
primarily on: 1) the systematic errors on the antiproton production and inelastic scattering 
cross sections; 2) the large uncertainty on the propagation parameters in the Galaxy and in 
the Heliosphere. 

Since several papers have been already published on this subject, we briefly mention 
here the main novelties of our approach. 

• We use, for the first time, all the relevant GR measurements taken during the same 
period and from the same experiment (PAMELA). 

• We use a recent computation [14] of the antiproton production and inelastic cross 
sections based on the newly available NA49 experimental data [15]. 

• We account for charge dependent solar modulation, including charge sign dependent 
drifts, by means of the recently developed Helioprop code [16]. 

Our approach goes beyond providing an updated range of allowed propagation setups to 
be used to propagate dark matter annihilation products. The widely used MIN-MED-MAX 
propagation models are defined as those models giving the maximal, median, and minimal 
supersymmetric (primary) antiproton flux and are compatible with B/G analysis [17]. Such 
set of models, however, cannot be confused, as often done in the literature (see however 
[18, 19]), with those providing the allowed range of secondary antiprotons. Indeed, while 
the primary antiproton flux is strongly dependent on the diffusive halo scale height and 
on the physical conditions in the Galactic center region, once the propagation parameters 
are set to fit the B/G, the secondary flux is almost independent on those quantities. Since 
the models providing the maximal/minimal secondary antiproton flux are dependent on the 
antiproton energy, we obtain for each energy bin the allowed flux compatible with GR nuclear 
measurements and the cross section uncertainties. We then constrain DM models using the 
minimum secondary antiproton flux evaluated in that way as background. 

Our results allow us to determine up-to-date conservative constraints on a representative 
sample of DM models and to compare these limits with the recent claim of a gamma-ray excess 
in the GG associated with DM annihilations in that region. 

In the last section, the preliminary results recently released by AMS-02 up to 450 GeV/n 
[20] are considered as well. In particular, we investigate the presence of an excess in the p/p 
spectrum at energies higher than previously measured by PAMELA. Given that the AMS-02 
B/G data are still preliminary, we do not aim at providing a dedicated statistical analysis. 
On the other hand, we select a single model providing the best fit of AMS-02 B/G data, 
and we tune the injection spectra to reproduce the proton and Helium data from the same 
experiment. Finally, we compare the predicted antiproton fluxes for this model with the new 
measurements. We discuss these results in the last Section. 
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2 CRs in the Galaxy 

2.1 CR propagation in the ISM 

We model CR propagation in the Galaxy using the publicly available numerical code DRAGON [21, 
22]. This code solves the diffusion-convection transport equation for all CR species in the 
Galaxy taking into account all the relevant processes, including energy losses and spallation 
due to the interaction with the ISM. We adopt the nuclear cross-section database taken from 
the public version of GALPROP ^ [23]. 

We assume cylindrical symmetry in the Galaxy. The propagation region is a classical 
cylindrical box with radial extension R = 20 kpc and height (L) varying in the range (2 — 16) 
kpc. The lower limit of the halo size is determined by the observed diffuse synchrotron emis¬ 
sion and by comparison between the computed low-energy secondary positron spectrum and 
PAMELA data [24, 25]. The upper limit is chosen in agreement with the available measure¬ 
ments of the ^^Be/®Be ratio. In the next section we will show that secondary antiprotons 
from CRs are almost unchanged by varying L in the given range, making the choice of the 
allowed domain not crucial in this context. We optimize the spatial resolution, as well as the 
time steps involved in the numerical algorithm, in order to evaluate the against PAMELA 
data with a precision <1%. 

We adopt a two-dimensional (r and z) propagation configuration since it is commonly 
used to study secondary CR antiproton production and also to keep the computational time 
short enough to scan over a large number of models. A more realistic three-dimensional 
approach could give a more accurate prediction of the CR fluxes, by taking the spiral arm 
distribution of CR sources and of the ISM gas into account [26, 27]. However we confirmed 
that these differences are much smaller than the other uncertainties under investigation in 
this work. 

We adopt the following standard expression for the diffusion coefficient D: 

Dip) = Do 13^ , ( 2 . 1 ) 

where p = p/Z is the rigidity of the nucleus characterized by charge Z and momentum p, and 
r] parametrerizes the low-rigidity behavior of D. While kinetic quasilinear theory predicts 
rj = 1 for the dependence of the diffusion coefficient on the particle speed, several effects may 
give rise to a different effective behavior (see e.g. [28, 29]). Bearing these considerations in 
mind, we leave r/ as a free parameter to be fitted against low-energy CR data. 

Diffusive reacceleration is parameterized in terms of the Alfven velocity va and we as¬ 
sume that this process takes place in the entire propagation region. In fact, this circumstance 
is expected if CRs are responsible for generating their own turbulence, e.g., by CR-induced 
streaming instability [30, 31]. Our approach is then different from semi-analytical methods 
where reacceleration is active only in the Galactic disk with height ~ 100 pc (see, e.g., [32]). 

We also account for a convective wind with velocity Vciz) vanishing at z = 0 and 
growing linearly with z with the gradient dVc/dz. Although dVc/dz can be as large as 
~ 100 km s“^ kpc“^ in the inner few kpcs of the Galactic disk [33, 34], the observables 
we are considering in the present analysis do not probe that region. On the other hand, 
the local convective velocity is severely constrained by ^^Be/®Be measurements [35] and, for 
this reason, we require dVc/dz to be less then 10 km s“^ kpc“^. By making this choice 

^http://galprop.Stanford.edu 
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we introduce another difference with respect to semi-analytical approaches where a constant 
convection velocity is adopted (see e.g. [32]). 

For the CR astrophysical source distribution we consider the following form: 

Qi{Ek,r,z) = fsNR{r,z) qo,i ( , (2.2) 

\ Po J 

Our spatial source distribution (/snr) follows the Galactic supernova remnants profile 
inferred from pulsars and stellar catalogues as given in [36] . We allow different source spectral 
indexes for the different nuclear species {'ji), and we fix their values - together with the species 
relative abundances ((?o,i) “ against recent experimental data. 

Proton and helium spectra exhibit a change of slope at a rigidity between ~ 200 and 
300 GV [11]: We model it by assuming a spectral break in the injection slopes of these nuclei. 
These features have relevant implications for the secondary antiproton flux. An alternative 
possibility for the origin of this break is a change in the ISM turbulence power spectrum 
and hence in the diffusion coefficient. In this context, high-energy GR antiprotons could be 
successfully used to discriminate among the two interpretations [13], but in the energy range 
where antiprotons are currently measured the two scenarios are quantitatively equivalent. 


2.2 Secondary antiproton production in the ISM 

The source term of secondary antiprotons produced from the interaction of CR nuclei with 
the ISM is given by the convolution of the antiproton production differential cross-section, 
da/dE^, and the interstellar CR nuclei differential energy flux, d^/dE^. 

The differential p production rate per volume and energy takes the form: 


Qp{Eu)= E E 

2=H,He j=H,He 


poo 

/ dE', 




(2.3) 


where E'^^ and E^ are the kinetic energies per nucleon of the incoming nucleus (with 
threshold energy E^h = 6mp) and the outgoing antiproton, respectively, n* denotes the 
interstellar gas density. 

Existing parameterizations of the antiproton production cross-sections [37] are mainly 
based on experimental data earlier than 1980. Recently, two works ([14, 38]) have improved 
these old computations making use of new precision data from the NA49 experiment at 
CERN [15]. In particular, in [14] the authors extended previous calculations by including the 
new laboratory measurements and a careful treatment of antiprotons arising from antineutron 
and hyperon decay. By using a revised approach to proton-helium scattering, the authors of 
[14] are also able to compute the production cross-sections for the processes: p + p ^ pX, 
p + He ^ pX, He + p ^ pX e He + He —)• pX, where the first particle is the impinging 
primary CR while the second one is the target interstellar material. 

We make use of their results for the antiproton cross sections, as well as their estimate 
of the related errors, in order to evaluate the nuclear uncertainties on our prediction of the 
secondary antiproton fluxes in the Galaxy. Finally, we also checked that our conclusion are 
unchanged by adopting the parameterization for the p + p scattering proposed by [38] and 
based on the full set of available measurements. 
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2.3 Antiprotons from dark matter annihilations 

Other than spallation of CRs on the ISM nuclei, antiprotons can be produced in the Galaxy by 
DM in pair annihilation or decay events. Antiproton measurements provide an invaluable tool 
to constrain such primary contribution, since the ratio between DM signal and background 
from standard astrophysical sources is usually much larger in the antiproton channel with 
respect to all other indirect detection methods. 

We compute the DM contribution to the antiproton flux as described in several previous 
works, e.g., [13, 39]. In brief, we assume that the source function (Qdm) for the WIMP DM 
component scales with the DM particle number density times the probability of annihilation 
and the antiproton yield per annihilation: 

Q-DuiEk^r^z) = ^ (2.4) 

where {av) is the thermally averaged annihilation cross section and /ODM(ic) is the DM density 
profile as function of the galactocentric distance x = y/r'^ + z'^. 

The DM profile is only poorly constrained by direct observations and its shape is usually 
inferred from N-body simulations of gravitational clustering. In our analysis, we adopt two 
spherically symmetric profiles: a standard Navarro-Frenk-White (NFW) [40], and a gener¬ 
alized NFW (gNFW) as defined, e.g., in [39]. The free parameters for the NFW profile are 
chosen following the analysis in [41], while for the gNFW profile we adopt a contracted profile 
(7 = 1.2) as suggested by the recent claim of a DM associated excess in the GC region [57]. 

For the choice of the annihilation channels, we focus here on two sample cases which have 
been recently investigated in connection to hints of DM signals (both direct and indirect), 
but potentially giving a sizable antiproton flux as well [13]. In particular, we consider as 
primary annihilation channels: DM DM —)• bb and DM DM —)• W^W~, and we take the 
corresponding antiproton yields from the PPPC4DMID [43, 44]. 


2.4 Propagation in the Heliosphere 

Low-energy (< 10 GeV) charged GRs are influenced by the solar magnetic field in the final 
stage of their propagation process. The main effect of the interaction with the heliosphere is 
a momentum decrease and the consequent alteration of the interstellar spectrum. 

This process is usually described in the context of the “force field approximation” [45]. 
According to this model, the energy spectrum d^T'OA/dEk^TOA of particles reaching the top 
of the atmosphere (TOA) with kinetic energy Fifc,TOA and momentum ptoa is related to the 
local interstellar spectrum (LIS) d^i,is/dEk as it follows 


d^TOA _ Ptoa 
dEk ,toa pIis dEk^his' 


Ek,Lis = Ek, toa + \Ze\(j), 


(2.5) 


where (p is the Fisk potential and parameterizes the kinetic energy losses. Within this ap¬ 
proach the effects of modulation for protons and antiprotons are identical. 

More realistic drift models, however, predict a clear charge sign dependence for the helio¬ 
spheric modulation of positively and negatively charged particles. Moreover, the modulation 
effect depends on the polarity of the solar magnetic field (SMF), which changes periodically 
every ~ 11 years. The SMF is observed with opposite polarities in the northern and southern 
hemispheres. At the interface between opposite polarity regions a heliospheric current sheet 
(HCS) is formed. The phenomenologically parameter a, known as the “tilt angle”, sizes the 
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Figure 1. The good models as discussed in section 3.1 in the S — Dq (left) and va — V (right) planes. 
Blue points correspond to models compatible with the PAMELA B/C ratio. Green points show the 
good models obtained using the PAMELA primary proton in addition. 


angular extension of the HCS oscillations. The magnitude of a depends on the solar activity 
and has a large influence on the intensity of the modulation. 

A more accurate description then may be achieved by means of dedicated numerical 
simulations of the CR transport in realistic models for the interplanetary magnetic field. 

In order to assess the impact of realistic model for solar modulation, we make use of 
the recently developed HelioProp code [16]. This code solves the equation describing CR 
transport in the heliosphere by means of the stochastic differential equation method [46, 47]. 
To characterize the model for solar propagation we have to specify the solar magnetic field 
geometry, the properties of diffusion, and those of winds and drifts. 

For the interplanetary diffusion tensor we assume H oc diag(A||, Ay^g), where the 
parallel and perpendicular components are set with respect to the direction of the local 
magnetic field. The CR mean free path parallel to the magnetic field as function of the particle 
rigidity p is parameterized as Ay = Ao(/9/lGeV)^(R/Rsun)~^) where B is the magnetic field 
and Rsun is its normalization value at the Earth position, for which we adopt Rsun = 5 nT. 
For the perpendicular diffusion we assume Ay^^ = Ay^g = 0.02A||, as results from numerical 
simulations [48]. 

Finally, for modeling the magnetic field geometry, winds and drifts associated to the 
antisymmetric components of the diffusion tensor we follow [16]. 

3 Method and results 

3.1 Selection of CR propagation models 

We focus our attention on the free parameters of our propagation model (see section 2.1) 
within the range given in table 1. We notice that the observables we are considering in our 
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Figure 2. The envelope of the secondary antiproton spectra computed with the different propagation 
models found to reproduce the B/C and primary spectra. 
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Figure 3. Left plot: The secondary antiproton flux computed for different values of the halo height 
L. Right plot: The normalization of the diffusion coefficient required, for each value of L, to reproduce 
the B/C ratio (the green line is used to guide the eye). 


analysis are not sensitive to L, we then use a fixed value L = 2 kpc and we evaluate the 
impact of different values for this parameter afterwards. 

For each model selected randomly in the parameter space described above, we fix the 
injection slopes and the source abundances for primary nuclei heavier than carbon by fitting 
CREAM data above 10 GeV/n [49] . For the carbon, helium and proton parameters, including 
the Fisk potential (/>, we fit the recent PAMELA Carbon [12] and proton [11] measurements 
at energy below break at ~ 200 GV. We assume Boron is entirely secondary. It was shown 
by [50] that, neglecting the production and acceleration of secondary nuclei inside SNRs, 
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Parameter 

Min value 

Max value 

Units 

V 

-1 

1 


Do 

0.1 

10 

10 ^^ cm^/s 

6 

0.2 

0.8 


dVc/dz 

0 

10 

km/s/kpc 

Va 

0 

100 

km/s 


Table 1. The free parameters of our propagation model with their associated range. The parameter 
L is taken fixed to 2 kpc. 


Observable 

Data points D 

# of parameters V 

XScr 

B/C 

18 

5 

29.79 

Proton flux 

71 

7 

98.21 

Carbon flux 

16 

7 

21.58 


Table 2. values corresponding to 3cr for the different observables. 


the S may be underestimated by a factor of ~ 5 — 15% (see also [51]). We checked that 
the Fisk potential gives an accurate description of modulated spectra compared against the 
more realistic predictions provided by the Helioprop simulations. Using a charge-dependent 
formalism for the modulation is relevant only when we compare differently charged particle 
spectra. We discuss this in detail in section 3.4. 

With the given set of diffusion and source parameters we are now able to calculate the 
B/C ratio. 

We identify a model as a good one, if it reproduces the B/C data as well as proton and 
carbon data within the Sa limits. 

In particular, we compute the x^’s against B/C (Xbc)’ proton (Xp) and carbon (xc) 
measurements for each propagation model and we accept the model if each x^ is smaller than 
the corresponding threshold (x^a) reported in table 2. The Xsa values reported in the table 
are calculated according to the distribution function with a number of degrees of freedom 
T = V — V with T) equals to the number of data points and V to the number of parameters. 
The 3 it threshold implies a 99.73% probability to get a smaller x^ value or, equivalently, a 
p-value of 0.27%. 

A similar strategy (using only Xsc) been adopted in [32] to select the models in 
order to determine the well-known MIN, MED and MAX setups. 

We repeat the procedure until N = 10^ good models are selected. 

In figure 1 we show where the selected models are located in the 5-Dq and VA-f] planes. 
We find that PAMELA data allow 5 to vary between 0.2 and 0.8 and this parameter is 
strongly anti-correlated with Dq. This anti-correlation can be explained as follows: An 
increase of 5 corresponds to a lower secondary-to-primary ratio at high energies, therefore 
Dq must be smaller in order to enhance the high-energy secondary production. The low- 
energy parameter 77 is strongly degenerate with the solar modulation potential and, as a 
consequence, the available data do not constrain 77 within the chosen range. By contrast, the 
Alfven speed parameter {va) is found to be bounded by ~ 30 km/s if primary spectra are 
included in the analysis and no low-energy breaks in the injection slopes are admitted. 
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Efc 

[GeV] 

V 

Do 

units 

<5 

VA 

[km/s] 

dVc/dz 

[km/s/kpc] 

Ip 

qc 

[xl03] 

ic 

0 

[GV] 

Min models 

1 

0.30 

3.32 

0.30 

32.2 

0.04 

2.58 

2.74 

2.53 

0.77 

10 

0.68 

2.85 

0.38 

28.6 

0.03 

2.54 

2.83 

2.48 

0.86 

100 

-0.16 

1.17 

0.75 

9.31 

6.78 

2.38 

3.40 

2.17 

0.57 

Max models 

1 

0.84 

0.85 

0.74 

0.52 

5.65 

2.40 

3.80 

2.18 

0.53 

10 

-0.92 

0.83 

0.68 

7.71 

4.05 

2.44 

4.03 

2.22 

0.54 

100 

0.60 

2.85 

0.23 

27.4 

6.88 

2.62 

2.95 

2.59 

0.75 


Table 3. Model parameters giving the minimum (maximum) contribution of secondary anti-protons 
at energy E = 1, 10, 100 GeV. 


The propagation model giving the best fit of the PAMELA B/C and proton data is 
characterized by the following parameters: Dq = 1-6, 5 = 0.41, va = 8.5, dVc/dz = 1.6, 
7 c = 2.56, 7 // = 2.47. 

3.2 Extreme models determination 

In order to evaluate the propagation uncertainty in the determination of the secondary an¬ 
tiproton flux, we show in hgure 2 the envelope of the secondary antiproton spectra computed 
with the propagation models selected beforehand. At lower energies the uncertainty band 
widens since more parameters are necessary to model CR propagation, while at larger ener¬ 
gies the main uncertainty comes from the poor determination of the 5 parameter. Our result 
is in agreement with [19]. 

We notice that the propagation model giving the best minimum (maximum) flux of 
secondary antiprotons is not univocally determined over the entire energy range 0.1—100 GeV. 
In fact, at the different energies the minimum (maximum) flux is associated with a different 
propagation setup. 

To give an example of the different models selected by varying the energy at which we 
evaluate the extreme fluxes, we provide in table 3 the model parameters associated with the 
minimum and maximum antiproton flux at 1, 10 and 100 GeV. 

Some trends emerge from this comparison. The minimal models selected at higher 
energies feature a larger S, while high reacceleration reduces the antiproton flux at low ener¬ 
gies. As expected the maximal models show the opposite behavior, indeed the hardest value 
allowed for 6 gives the maximum contribution to the antiproton flux at higher energies. 

It is important to remark here that our poor knowledge of the halo size does not affect 
these conclusions. To show this, we select the propagation model giving the B/C best ht 
and we test different values for L up to 16 kpc. In order not to lose the perfect agreement 
with the secondary over primary data, we increase the Dq value accordingly (see the right 
plot in hgure 3). As shown in hgure 3, different choices for L in this range do not affect our 
predictions for the secondary antiproton hux. 

Although in this paper we assumed a uniform value of <5 in the whole Galaxy, it was 
recently shown that diffuse 7 -ray data favor a scenario characterized by radially-dependent 
CR transport properties [52, 53]. In order to investigate the possible impact of that scenario 
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Figure 4. Comparison between propagation and nuclear uncertainties. Yellow band: Error on the p 
flux due to the uncertainty in the propagation parameters. Blue lines: The relative difference between 
the p flux computed using the fiducial cross section from [14] and: (dot-dashed) the maximal model 
from [14], (solid) the minimal model from [14], (dashed) the parameterization from [23, 37]. 


on our results, we computed the local secondary antiproton spectrum for the KRA.^ model 
considered in those papers finding a negligible correction. 

3.3 Antiproton production cross-section uncertainties 

We compare here the propagation uncertainties derived in the previous sections with those 
associated with the antiproton production processes. 

In figure 4, we show the relative difference between the minimum (maximum) secondary 
antiproton flux and that obtained using the best-fit propagation model. The corresponding 
region represents the uncertainty on the secondary flux associated with galactic propagation. 

We compare this uncertainty band with the relative differences associated with produc¬ 
tion cross sections. To this end, we compute secondary antiprotons with the new prescriptions 
recently proposed by [14] and we evaluate them against the traditional fitting relations given 
in [23, 37]. 

We find that nuclear uncertainties can be as large as 50% even at ~ 100 GeV, and 
are much larger below few GeVs. However, with the available GR data, the propagation 
uncertainties dominate over the entire energy range as shown in figure 4. 

Upcoming measurements (in particular, from AMS-02 [1], GALET [54], andlSS-GREAM [49]) 
are expected to significantly improve our knowledge of propagation parameters and then to 
reduce the associated uncertainties. In that situation, antiproton production cross sections 
will prevent us to provide predictions for the astrophysical backgrounds as accurate as the 
forecasted sensitivities. 

3.4 The role of charge-dependent solar modulation 

As pointed out in section 2.4, charge-dependent solar modulation can be relevant when the 
TOA antiproton flux is evaluated. Therefore, we compare here our predictions of the extreme 
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Figure 5. The envelope of the secondary antiproton spectra computed with the charge-dependent 
modulation {black lines) and compared with that one obtained with the force-field approximation 
{yellow band). 


fluxes based on the force-field approximation with those obtained with a charge-dependent 
modulation model. 

To modulate the antiproton flux in the charge-dependent scenario, we develop the fol¬ 
lowing strategy: 

• For each propagation model, we consider as free parameters: 1) solar magnetic field 
polarity; 2) a (HCS tilt angle); 3) Aq (normalization of the parallel mean free path); 4) 
5 (power-law slope of the heliosphere diffusion coefficient as function of rigidity). Solar 
polarity and a are fixed by the data-taking period, since they can be obtained by direct 
measurements [55] , while we determine Aq and b by fitting the predicted TOA proton 
flux against the low-energy PAMELA measurements. 

• We use the same set of parameters obtained from protons to modulate the LIS antipro¬ 
ton flux. 

In figure 5 we show the extreme antiproton fluxes as obtained with our charge-dependent 
modulation model. We immediately notice that the more detailed treatment of solar modula¬ 
tion does not impact significantly on the determination of the minimum flux, while differences 
up to ~ 20 — 30% are shown for the maximum one (see also [39] for a more detailed discussion 
on this issue). 

To understand this result, we point out that the models giving the minimal antiproton 
TOA fluxes are the ones with the largest LIS proton spectrum and, therefore, they need a 
stronger modulation to reproduce the data. Stronger modulation is associated with a larger 
Fisk potential or a smaller heliospheric parallel mean free path. In this situation, diffusion 
dominates over charge-dependent drifts and the modulation of particles with different charges 
becomes equivalent. 
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Figure 6. Uncertainties on the flux of primary antiprotons originating from DM annihilation due to 
the propagation parameters (yellow band) and to the halo size (red line). For the latter, the ratio 
between the two extreme cases L = 16 kpc and L = 2 kpc is considered 


4 Conservative limit on DM models from antiproton data 

In this section we derive the most conservative - with respect to all the background un¬ 
certainties discussed before - constraints on the DM annihilation cross section for the DM 
WIMP scenarios introduced in section 2.3. 

To propagate DM antiprotons, we choose L = 2 kpc since it is the minimum value 
compatible with synchrotron diffuse emission observations [24, 56]. We note here that, while 
the actual value of L is irrelevant for the secondary antiprotons (see Fig. 3), DM antiprotons 
can change significantly and, in fact, this parameter is the most important one to evaluate 
this contribution. In particular, thinner halos underproduce the DM p flux, and therefore 
L = 2 kpc corresponds to the minimum flux expected from a given DM model (see [13] for 
a more detailed discussion) and, for that reason, to the less stringent bound. 

In order to quantify these statements, in Fig. 6 we show the ratio between the antiproton 
flux corresponding to L = 16 kpc and L = 2 kpc, compared to the uncertainty band due to 
the poor knowledge of the propagation parameters. The plot clearly shows how, as far as 
primary antiprotons from DM annihilations are concerned, the uncertainty on the halo size 
is strongly dominant. 

We also remark that we can safely neglect the charge-dependent effects in the determi¬ 
nation of the minimum background (see section 3.4). 

In order to get the most conservative bound with respect to the propagation setup, 
a naive strategy would be the following: Starting from a minimal background, the DM 
component is added until a best ht to the data is reached, and then the DM cross section 
is increased until the quality of the fit is worsened up to the 2 ct level. However, although 
the models providing the minimal background are compatible with the B/C within 3a, they 
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Figure 7. Antiproton bounds on DM annihilation rate. Red lines: bb channel for NFW profile for 
different assumption for the secondary p production. Blue lines: the same for the WW channel. The 
results obtained with a gNFW profile with 7 = 1.2 are indistinguishable from the NFW ones. 


do not necessarily provide a satisfactory fit of the antiproton flux. The addition of a DM 
component could still leave unexplained the antiproton flux at energies above the DM mass. 

Hence, we determine the 2 ct exclusion contour in the plane (ut-dM) as follows. 

For each good propagation model (based on B/C and nuclear data), we first compute 
the secondary background and the DM flux for a given DM mass, channel and profile. We 
then find the best-fit value for the annihilation cross section. 

If the combination of background and DM (computed assuming the best-fit cross section) 
satisfies the following condition: 

x" {(Thest&tv) < xio5 (-F = 23 - 1 = 22) ~ 33.92 (4.1) 

(where the number of degrees of freedom iF is computed taking into account 23 data 
points and 1 parameter), then we retain the model. In other words, we reject propagation 
models with a p-value lower than 5%. 

If the model passes this test, we compute the value of {av) above which the fit worsens 
beyond the 2a level with respect to the best fit. More precisely, we set the following threshold, 
corresponding again to a p-value equal to 5%: 

X^ (o-max^^) - X^ (o-bestfit^^) < Xao5 (-F = 1) ~ 3.84 (4.2) 

where now the only degree of freedom is the annihilation cross section. 

Finally, a scan over the models is performed: The limit value shown in Fig. 7 is therefore 
obtained by taking the maximum cross section (with respect to the different propagation 
models) satisfying the bound described above for given DM mass. 

We remark that this strategy allows to obtain the bound by propagating secondary and 
DM antiprotons consistently with the same propagation model. 
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In the same plot, we compare the limits we obtain when the background is computed 
with the hducial model for the p production cross-section and with their minimal and max¬ 
imal values. Because of the strategy adopted to derive them, the case in which the most 
conservative limits are obtained with the minimal production cross sections is not always 
prevalent: Indeed, it may happen that larger production cross sections select a different 
subset of propagation models to reproduce the p measurements in combination with DM. 
Some of the propagation models allowed with fiducial or maximal cross-section, but not with 
the minimal case, may produce a lower DM p flux leading to a larger 2a annihilation cross 
section. 

In figure 7 the reader can see our results for the maximum allowed annihilation cross 
section for the bb and W~^W~ annihilation channels. The maximum allowed cross section we 
find for bh is around one order of magnitude larger than what the authors of [39] found for 
the charge-symmetric modulation case. The main difference in our approach is to make use 
of a broader analysis of the propagation and nuclear uncertainties in order to determine the 
background. 

Our results can be now compared with the DM interpretation of the recently claimed 
signal in the gamma-ray channel located in the inner few degrees around the GC [42]. 

In [42] the authors show that a DM particle with mass ~ 43 GeV annihilating into bh 
with a cross section (av) ~ 2.2-10“^® cm^s“^ (for the Inner Galaxy analysis) and distributed 
according to a gNFW prohle with 7 = 1.2 can accomodate the anomalous excess. 

The detailed analysis reported in [57] provided a better quantification of the system¬ 
atic uncertainties affecting the proposed signal; more recently, the Fermi-LAT collaboration 
released preliminary estimates of the energy spectrum of this excess, based on four qualita¬ 
tively different background models [58] . A wide set of DM masses and annihilation channels 
are compatible with these new analyses (see, e.g., [59]). It has been also shown that, in the 
context of the Minimal Supersymmetric Model framework, these candidates are not in ten¬ 
sion with LHG or direct detection constraints. In figure 7 we compare our hndings with the 
favored regions of annihilation cross sections connected to the GC excess as reported in [57] . 

The bottom line of this analysis can be summarized as it follows: Although some fiducial 
models of CR propagation would produce strong tension with the DM interpretation of the 
GC excess (see, e.g., [39, 60]), given the large uncertainties on the propagation parameters 
(for the secondary p) and on the halo height (for the DM p), the antiproton channel cannot 
be invoked to conclusively exclude this hypothesis. 

5 Discussion on AMS-02 data 

In this section we focus on the recently released AMS-02 data, including protons [61], and 
preliminary helium, B/C and p/p ratio [62], with energy range extending to 450 GeV. 

In particular, we take a closer look at the new impressively accurate data on the p/p 
ratio and we attempt to evaluate their compatibility with the other hadronic observables. 
Given the preliminary nature of the released data we do not attempt a statistical analysis 
of the uncertainties associated with propagation. In this perspective, the final release of the 
secondary/primary measurements, when systematic and statistical errors are fully accounted 
for, will be crucial. 

A propagation model chosen among those considered in section 3.1, and compatible 
with preliminary B/G measurements, is shown in figure 9. The propagation parameters are: 
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Figure 8 . Our reference model is compared to AMS-02 proton [61] {left plot) and helium {right plot) 
data. With the dotted and dashed lines we show the minimal and maximal breaks compatible with 
the hardening measured by AMS-02. 
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Figure 9. Our reference model compared to AMS-02 preliminary B/C data. Solid line: the TOA 
spectrum modulated with (f = 0.6 GV; dotted line: the LIS spectrum. 


Do = 1.5, 6 = 0.42, VA = 27, dVc/dz = 14, 7 c = 2.46, 7 // = 2.44/2.31. For comparison, the 
same value for S was found by [51] using the same datasets. 

Remarkably, the predicted B/C ratio reproduces the AMS-02 data over more than three 
orders of magnitude in energy. It is worth noting here that the 6 required by the new high- 
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Figure 10. Our reference model compared to AMS-02 preliminary p/p data. Blue solid line: the 
p/p spectrum computed with the fiducial cross sections from [14], with the optimal hardening in the 
proton and helium injection spectra. Dotted and dashed lines: the p/p spectrum computed with the 
minimal and maximal hardening in the proton spectrum as in Fig. 8. The blue band reports the 
uncertainty associated to the production cross sections. 


energy measurements is in perfect agreement with the best-fit value obtained in our earlier 
statistical analysis [ 6 ] , based on the available high-energy measurements preceding PAMELA 
and AMS-02 releases. 

We also tune the proton and helium injection slopes to accomodate the AMS-02 data. 
For the protons, we also consider the minimal and maximal injection slopes at high energy 
compatible with the data. The reader can see the comparison with the new datasets in 
figure 8 . 

Armed with a model fully consistent with all the preliminary nuclear observables, we 
can finally compare our prediction for the p/p ratio with the data. 

In figure 10 we show this comparison. The computation of the secondary flux is per¬ 
formed using the fiducial value of the cross sections provided by [14], and the associated 
uncertainty is shown as a blue band. 

We conclude that, even without considering all the relevant uncertainties associated 
with propagation or injection slopes, our predictions for the p/p are in good agreement with 
the preliminary data in the entire energy range. Our findings are then in agreement with the 
conclusions of [63], although our analysis relies on the B/C data from the same experiment 
for the assessment of the propagation model. 

6 Conclusions 

We presented a revisited study of the dominant uncertainties in the determination of the CR 
secondary antiproton spectrum. 
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By performing a scan over the parameter space relevant for CR propagation, we iden¬ 
tified a set of models compatible with B/C, proton, helium and carbon data provided by 
the PAMELA experiment. We were then able to bracket the minimum and maximum sec¬ 
ondary antiproton fluxes constrained by local observables and we compared the associated 
uncertainty band with the errors related to the production cross sections. It is the first time 
that such analysis has been performed by using comprehensive numerical simulations of CR 
propagation in the Galaxy and the Heliosphere. More importantly, we used for the first time 
a complete set of measurements from the same experiment: Using consistent data from the 
same data-taking period allowed us to reduce the uncertainties due to solar modulation. 

Similarly to previous results, we found that the determination of the (almost unknown) 
diffusion halo height is irrelevant for the computation of the secondary antiproton flux since 
this parameter is degenerate with the diffusion coefficient normalization Dq. In addition, 
we found that using the recent PAMELA data, the uncertainty on the propagation model 
dominates over the nuclear ones. 

Our result has important implications for the indirect search of primary p from DM 
annihilations in the galactic halo. Therefore, we provided the most conservative - with 
respect to the mentioned effects - constraints on the annihilation rate for some popular DM 
models recently investigated in connection to hints of DM signals in other detection channels. 

Our method may be taken as a reference procedure to be exploited when the final 
measurements for all the relevant channels are published by the AMS-02 collaboration. 

At the moment, the preliminary release by the AMS-02 collaboration of nuclear data 
does not permit to perform a statistical analysis. Nevertheless, we found that the model in 
agreement with AMS-02 proton, helium, and B/C data is compatible with the p/p spectrum. 
Therefore, we do not report any significant anomaly in this observable. Our result is then 
consistent with the conclusions presented in [63]. 
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